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ABSTRACT 

In a thermally bistable medium, cold, dense gas is separated from warm, 
rareified gas by thin phase transition layers, or fronts, in which heating, radiative 
cooling, thermal conduction, and convection of material are balanced. We calcu- 
late the steady-state structure of such fronts in the presence of magnetic fields, 
including the processes of ion-neutral drift and ion-neutral frictional heating. We 
find that ambipolar diffusion efficiently transports the magnetic field across the 
fronts, leading to a flat magnetic field strength profile. The thermal profiles of 
such fronts are not significantly different from those of unmagnetized fronts. The 
near uniformity of the magnetic field strength across a front is consistent with 
the flat field strength-gas density relation that is observed in diffuse interstellar 
gas. 

Subject headings: diffusion — ISM: structure — methods: numerical — MHD 



1. INTRODUCTION 

The low- and intermediate-temperature parts of the interstellar medium (ISM) consti- 
tute a thermally bistable medium that results from the balance between radiative heating 



and cooling as well as heating by cosmic rays ( iField et al. Ill969[ ) and photoelectric heating 



from PAHs (jWolfire et al. Ill995l ). the dominant heating source. The two stable phases are 
referred to as the cold neutral medium (CNM), having T CNM ~ 10 1-2 K, and the warm 
neutral medium (WNM), with T WNM ~ 10 3 ~ 4 K. The de eree to which magnetic fields are 
frozen into this interstellar gas is parameterized by the magnetic Reynolds number, Reu, 
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and the ambipolar Reynolds number, Re ad (IZweibel fc Brandenberg 1119971 ). The magnetic 
Reynolds number is given by the ratio of the Ohmic diffusion time to the dynamical time, 
and for ISM parameters is of order 10 15 — 10 21 . The ambipolar Reynolds number, given by 
the ratio of the ion-neutral drift time to the dynamical time, is many orders of magnitude 
smaller and may approach unity in dense molecular gas. Based on these estimates, one would 
expect that magnetic fields should be well coupled to both the ionized part of the gas and 
to the neutrals for all but the most dense or low column density clouds. 

Under ideal magnetohydrodynamic (MHD) conditions, such as those indicated above, 
one might expect a strong correlation between magnetic field strength and density. If the 
relationship is expressed as B oc p x and we ignore diffusion, for flows directed transverse 
to the field we have x — 1> whereas for field aligned flows x = 0- The median magn etic 
field strength in the CNM has been measured at B ~ 6 pG ( IHeiles fc Troland 1 120051 ) . If 
the field was frozen in we might expect to detect much smaller field strengths in warmer, 
lower density gas, but instead it is found that the field strength in other ISM components is 
similar to that of the CNM. This was demonstrated by measurements of the Zeeman effect 
over the density range 0.1 c m" 3 < n < 100 cm" 3 tha t yielded a flat magnetic field strength- 



gas density (B — p) relation ( iTroland fc Heiles 1119861 ). The most obvious explanation for this 
relation is that motions are aligned with the magnetic field. However, this has been argued 
against in two ways. First, in order for a magnetic field to collimate a flow in this manner 
it must dominat e the turbulent energ y density, but the field strength is less than or equal 



to equipartition ( IHeitsch et al. 1 120041 ) . Second, the accumulation length for the formation 



of giant molecular clouds is of order a kiloparsec and may be too large a scale over which 
to expect coherent flows (IMestel Ill985l ). Thus, the flat B — p relation may be indicative of 
magnetic diffusion. 

Among the mechanisms that have been proposed to account for the flat B — p relation 
in the diffuse ISM are turbulent ambipolar diffusion (Zweibel 2002; Heitsch et al. 2004), 
decorrelation due to MHD waves fjPassot & Vazquez-Semadeni 1120031 ). and turbulent mag- 
netic reconnection (e.g. Santos-Lima et al. 2010). These dynamical studies argued that 
ambipolar diffusion alone was not sufficiently fast to transport magnetic flux over the large 
scales under consideration and so invoked turbulence to enhance transport. However, a 1-D 
two-fluid dynamical study of the thermal instability as a formation mechanism for diffuse 
clouds showed ambipolar diffusion to efficie ntly transport magn etic field such that the ob- 
served B — p relation could be reproduced (jlnoue et al. 1 120071 ). The work presented here 
complements those findings but is also a significant departure from that and the other cited 
examples as we shall consider the actual transitions from one phase to another. Our ap- 
proach is advantageous in that we control the diffusive processes and are not hampered by 
numerical diffusion. The hydrodynamic structure of CNM/WNM transitions has already 
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been presented (Inoue et al. 2006; hereafter IIK06), but the effects of a magnetic field have 
not previously been studied. In the case of a magnetic field orthogonal to a transition layer, 
the field has no effect on the structure as it does not exert any force or modify thermal 
conductio n in the direction of th e temperature gradient, although it has a large effect on 



stability (IStone fc Zweibel II2009I ). In this work, we consider the case of a magnetic field 
tangential to a transition layer in a simple 1-D geometry and include ion-neutral drift as the 
magnetic diffusion mechanism. 

Our paper is organized as follows: in §2 we present our numerical method for calculating 
the structure of a phase transition layer for a given initial density and magnetic field strength. 
In §3 we discuss the effect of ambipolar drift heating on the two-phase structure of the neutral 
ISM. In §4 we show a selection of our ambipolar diffusion-mediated front solutions, and 
include a brief discussion of the flux-freezing approximation. In §5 we discuss the physical 
significance of our results, and in §6 we summarize our findings. 



2. METHOD 

We consider the scenario of a phase transition layer, or front, separating two uniform 
media of different densities and temperatures in a simple one- dimensional geometry with 
x as the direction of variation. A uniform magnetic field is tangential to the front such 
that B = B(x)z. The geometry is illustrated in Figure [TJ We assume a steady-state and 
ionization equilibrium, and work in the reference frame of the front. These assumptions 
shall be justified in §5. In order to calculate the physical structure of a front we consider 
six variables: pressure (p), density (p), bulk velocity (v), plasma velocity (v p ), magnetic 
field strength (B), and temperature (T), that are described by the following five equations, 
namely the equation of state: 

RpT 



the continuity equation: 
the momentum equation: 

the induction equation: 



dB d 



dt dx 



(v P B), (4) 
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and the energy equation: 



7 R dT dp d dT 

-P-TT ~ -77 = ( 5 ) 



7 — 1 p dt dt dx dx 

where 7 is the adiabatic index, R is the molar gas constant, p is the mean molecular weight, 
k is the thermal conductivity, pC is the cooling function (which includes ambipolar drift 
heating), and djdt = d/dt + vd/dx. In the approximation that the plasma and neutral 
fluids are well coupled, and the neutral density dominates, the plasma velocity may be 
written as the sum of the drift velocity, vd = Vi — v n , and center of m ass velocity, v « v n , 



such that v p ~ v + Vp, where the drift velocity is given by (IShu 1119831 ): 



vd = , (6) 

CPiPnlAD 

with the drag coefficient for collisions between ions and neut rals given by 7,4 p = < av >i n 



I (mi + m n ) cm 3 s 1 g 1 , where < av >i n = 2 x 10 9 cm 3 s 1 ( iDraine et al. 1119831 ) 



Assuming a steady-state (and having already dropped the y and z dimensions), inte- 
grating Equations ©, (EJ), and (JSJ) with respect to x yields the following conservation laws: 

3 = f>v, (7) 
B 2 

M B = pv 2 +p + — (8) 

87T 

7 R.dT dp d dT 

:3-r: - Vj- = tt-^tt- - pC, (9) 



7 — 1 /i dx dx dx dx 
where j is the mass flux and Mb is the total energy density. To solve Equation (J9]) we require 
expressions describing the evolution of the flow speed and magnetic field strength, including 
the process of ambipolar diffusion. An equation for the flow speed is obtained by taking the 
derivative of the total energy density, Equation OH]), to obtain: 



(10) 



dx 47rj(RT — pv 2 ) 

An equation describing the magnetic field strength is derived by substituting B = B(x)z 
into Faraday's law in 1-D and using Equation (@| to yield: 

dB z dE v 

~dT = ~ C ~dx ^ v = ~ X v = Vpx = constant - 

Substituting the plasma velocity, with the drift velocity given by Equation (jSJ), into Equation 
(ITTj) we obtain: 

B 2 dB 

vB - — = cE. 12 

^PiPnlAD dx 
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This is a first-order ODE with one parameter, cE. Mathematically, cE can take any value 
since it is a constant of integration. However, we will argue at the end of this section that 
physical considerations of the magnetic field strength and ambipolar heating across a front 
serve to greatly reduce the cE parameter space. 

Equations (Q, ({TO]) , and 012p . and the definition z = dT/dx yield a system of four ODEs 
for T, B, and v that apply for any functional form of conductivity and cooling function. In 
the gas states studied here, con ductivity is d ominated by neutral atoms such that k = 



2.5 x lO 3 ^/ 2 ergs s" 1 K" 1 cm 1 (IParker II1953T ). The cooling function is written in full as: 



p£ = n[nA-(r PA H + r AD )]. (13) 

We take the simple functional forms used by IIK06 for Lya and [C II] radiative cooling such 
that: 

91 / -118400 K\ 97 /-92 K\ , , 

A = 7.3 x 10~ 21 exp + 7.9 x 10" 27 exp ergs s^cm" 3 , 14 

p yT + 15001^ y T J ' v ; 

and for photoelectric heating: 

r PAH = 2 x 10~ 26 ergs s" 1 . (15) 
Heating by ion- neutral friction is represented by Tad'- 

2 1 ( B dB \ 2 -l -3 

nT AD = piPn^ADVo = - — i- ergs s cm (16) 

PiPnl AD y^TC dx J 

(Scalo 1977, Padoan et al. 2000), where the density of neutrals is given by p n = p n mHn, 
and the density of ions by p, = pim p n e . We compute the ionization fraction using: 

n e ( 4 s T°- 845 \ 

— = 1.19 x 10" 4 - 1.36 x 10" 8 + 

n H \ n H ) 

^0.845 jil.69 \ ' 

1.42 x 10~ 8 + 2.72 x 10~ 8 + 1.85 x 10~ 16 ^^ (17) 



flFerriere et al. Ill988h 



In solving this system we choose initial values for the density and magnetic field strength 
and impose the following boundary conditions: 

T(x = x x ) = T u T(x = x 2 ) = T 2 , (18) 
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where x\ and 0:2 represent the left- and right-hand boundaries, respectively, and T\ and T 2 
satisfy thermal equilibrium at these boundaries. T\ is found by solving Equation ffl3|) for the 
chosen initial value of the density at x±, and T 2 is the temperature obtained by integrating as 
far as X2, where the size of domain is chosen such that T 2 will also satisfy thermal equilibrium. 
For our third and fourth conditions, given by Equation (fl9|) . we impose zero temperature 
gradient at both boundaries. Finally, we set the value of the initial magnetic field strength 
gradient, \dB / dx\ Xl , as this controls the amount of ambipolar heating in a given front model. 
As we will show in §4, the choice of the initial field strength gradient affects the structure of 
the front. Given that we set five boundary values but have a system of only four ODEs, we 
thus set up an eigenvalue problem in which the mass flux, j, is the parameter to be adjusted 
to find a self-consistent solution. 

The numerical method we employ is that of shooting, in which the integration is per- 
formed with an initial guess for j, the resulting boundary values compared to the desired 
conditions, and j adjusted accordingly so that the integration can be repeated as necessary 
until the right-hand boundary conditions are satisfied to within some chosen tolerance. We 
find that the degree to which thermal equilibrium is satisfied at the right-hand boundary 
depends on the size of the domain, which should be adjusted to achieve optimum results. 
For cases in which ambipolar drift heating does not dominate it is possible to satisfy ther- 
mal equ ilibrium to better than one part in 10 5 . We use a 5th-order adaptive Runge-Kutta 



scheme (IPress et al. 1119921 ) and adjust the eigenvalue according to the secant method. When 
appropriate bounds are chosen our method converges to a solution quickly, requiring of order 
ten iterations. Note that we always integrate from the cold medium to the warmer one. 

We close this section with a brief discussion of the initial magnetic field strength gradient 
boundary condition and the parameter cE. In setting up our initial conditions, instead of 
choosing the value of cE directly we instead set the initial value of the magnetic field strength 
gradient, \dB/dx\ Xl . This implies the value of cE, which is kept constant across the domain, 
as we may evaluate it by substituting our initial conditions into Equation ffl2l) . Note that 
the value of cE will change with each iteration of the shooting method because it depends 
on the bulk velocity, which is adjusted according to the secant method. For all density and 
magnetic field strength initial conditions there is some minimum value of cE below which 
the magnetic field strength gradient is positive throughout the domain. The outcome for 
choosing an initial field gradient that yields a value of cE below this minimum would be a 
larger magnetic field strength in the warm medium than in the colder one. However, if one 
imagines an evaporating cool cloud with the assumption of frozen-in magnetic field lines, 
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this does not seem like a physically reasonable scenario as the field lines will become further 
apart as the cloud expands. Furthermore, if the value of cE is too large, ambipolar drift 
heating may dominate over photoelectric heating making it increasingly difficult to satisfy 
thermal equilibrium at the far boundary, implying that a front can no longer exist. We 
demonstrate quantitatively the effects of \dB/dx\ Xl , and hence cE, in §4, but do not refer to 
cE explicitly in the rest of the paper. 



3. EFFECTS OF AMBIPOLAR DRIFT HEATING ON TWO-PHASE 

STRUCTURE 

The two neutral phases of the ISM are enabled by the balance of radiative cooling 
and heating by, in this work, photoelectric heating and ion-neutral friction. We present 
the equilibrium state of the cooling function, p£(n, T) = 0, in Figure El with pC given by 
Equation ( II 3p . The solid line shows the case in which there is no ambipolar drift heating, for 
which IIK06 report that a two-phase structure is possible for 10 2 8 K cm -3 < p/ks < 10 41 
K cm -3 . The other lines illustrate the effects of increasing the ambipolar drift heating rate 
at a fixed magnetic field strength of 3 pG. Although we have already shown Tad to be a 
function of the density and field strength, for the purposes of this plot we have set it to be 
a constant fraction of the photoelectric heating rate, Tpah- Increasing the total heating rate 
serves to increase the pressure at which two phases can co-exist: the minimum pressure at 
which the cold phase can exist, and the maximum pressure at which the warmer phase can 
exist, both increase. In fact, the pressure range over which two phases can exist becomes 
larger as the total heating is increased. For example, for the r^D /Tpah = 0.50 case plotted 
in Figure [21 two-phase structure is possible for 10 30 K cm~ 3 < p/ks < 10 4 ' 3 K cm~ 3 , and for 
the Tad/Fpah = 1-00 case the pressure range is 10 31 K cm -3 < p/k B < 10 44 K cm -3 . We 
can understand the shift towards lower densities as follows: increasing the heating increases 
the temperature, so it must decrease the density. The upshift of the equilibrium to higher 
pressures also reflects the increased heating. 



4. FRONT SOLUTIONS 

The characteristics of a front are determined by its thermal pressure. There exists a 
"saturation pressure" at which heating and cooling are balanced within a front (Zel'dovich & 
Pikel'ner 1969, Penston & Brown 1970). If A and T can be written as functions of pressure 
and temperature (where T is the total heating rate), this pressure may be calculated by 
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solving the integral (llnoue et al. 1120061 ): 




rai(nA - T)dT = 



(20) 



and substituting for n using the equation of state, Equation (0Q). For the hydrodynamic case 
{T ad — 0) IIK06 obtain p sa t/kB = 2612 K cm -3 (where ks is the Boltzmann constant), 
which, by solving Equation 013 j) . implies an initial density of n = 106.08 cm -3 and hence an 
initial temperature of T = 24.63 K. If the thermal pressure exceeds this value of p sa t a fluid 
element passing through the front experiences net cooling, so we have a condensation front. 
If instead the thermal pressure is less than the saturation value a fluid element experiences 
net heating, so we have an evaporation front. 

In this section we demonstrate the effect of ambipolar diffusion on the saturation pres- 
sure and present our ambipolar diffusion-mediated front solutions. We also argue that the 
flux-freezing approximation is not accurate for steady-state thermal fronts. 

4.1. Effects of Ambipolar Drift Heating on Saturation Pressure 

The saturation pressure is altered in the presence of a magnetic field due to ambipolar 
drift heating. The integral given by Equation ( 120]) cannot be solved analytically when Tad is 
non-zero, so instead we use our shooting method, as discussed in §2, to find the initial density 
that yields a static solution as a function of the initial magnetic field strength gradient. The 
results for initial magnetic field strengths of 1, 3, and 5 /iG are shown in Figure El Note 
that the field gradients are actually negative, as we anticipate the magnetic field strength 
to decrease with increasing temperature, and we refer to the absolute magnitude of the 
quantity, which we give in units of /xG pc _1 . 

As \dB/dx\ Xl is increased the saturation density and pressure for all magnetic field 
strengths initially decrease until a sufficiently large value of \dB/dx\ Xl is reached, after which 
the density and pressure both increase. Therefore it is possible to have two different fronts 
at the same saturation pressure. This nonmonotonic behavior may be understood by solving 
Equation (120]) for the saturation pressure using the equation of state, Equation (JTJ), to obtain: 



This shows that increasing the total heating rate, T, tends to increase the saturation pressure. 
But increased heating also tends to drive up the temperature, which for CNM temperatures 
greatly increases the cooling rate, A, and according to Equation ( 12 ip this decreases the 



Psat 



k B 



rT 2 kA, t 



(21) 
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saturation pressure. For example, Figure [3] shows that if Bq = 5 /iG and \dB/dx\ Xl = 300 /iG 
pc _1 the CNM temperature is increased from 24 to 28 K. According to Equation (|14p . this 
results in a greater than 70% increase in the cooling rate, A. Such a large increase in cooling 
requires a lower density and a lower saturation pressure for equilibrium to be maintained. 
This effect dominates as long as the heating and cooling rates, T and A, are not too large 
and is the reason for the dip in the saturation pressure seen in Figure |3j An inflection point 
is not observed in the saturation pressure in the 1 /iG case, the reason being that at higher 
values of \dB/dx\ xl (and hence larger ambipolar heating rates) the magnetic field profile is so 
steep that a thermal equilibrium phase cannot be reached before the magnetic field strength 
becomes negative. In such instances the temperature on the cold side is still well within the 
range of CNM values, so it is not the medium being overheated that prohibits physical front 
solutions. 

We present example saturation fronts having a thermal pressure of pth/ks = 2500 K 
cm -3 and an initial field strength of 3 /iG, but with different initial values of \dB/dx\ Xl , 
in Figure HI Note that the values of \dB/dx\ Xl given represent the largest gradients at 
any point throughout the front. The magnetic field strength gradients of all the fronts we 
present quickly relax to become much smaller than the initial values that we impose. The 
front having the larger value of \dB/dx\ Xl has a higher ambipolar drift heating rate and 
connects a lower density, higher temperature CNM with a higher density, cooler WNM than 
the static front with the lower heating rate. The front with the lower ambipolar drift heating 
rate is the most diffusive, which is illustrated by its flatter magnetic field strength profile. 
This is also indicated by the ratio of the field strength to the number density, which shows 
a larger variation across the domain than the same quantity for the static front with the 
higher heating rate. 



4.2. Ambipolar Diffusion-Mediated Front Solutions 

As stated at the beginning of §4, in the hydrodynamic case a static front has a thermal 
pressure of Psat/ks = 2612 K cm -3 , which corresponds to an initial density and temperature 
of n = 106.08 cm -3 and T = 24.63 K, respectively. To demonstrate the effects of ambipolar 
diffusion we present several front models having this same initial density and an initial 
magnetic field strength of 3 /iG in Figure The initial temperature changes according to 
the ambipolar drift heating rate, so is not the same as in the hydrodynamic case. The 
different models correspond to various initial magnetic field strength gradients, \dB / dx\ Xl , 
where a larger gradient corresponds to increased heating. The properties of the phases 
connected by these fronts are listed in Table 1. 
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The overall shapes of the temperature profiles, shown in the top left panel of Figure 
are fairly similar with the main differences being the temperature gradients on small scales 
and the final temperatures of the warm phases becoming lower as \dB/dx\ Xl is increased. 
The main effect of increasing \dB/dx\ Xl is that the size of the integration domain required to 
reach thermal equilibrium at the right-hand boundary becomes smaller, due to the increased 
ambipolar drift heating. In fact, the lowest \dB/dx\ Xl profile shown here is very similar 
to the hydrodynamic solution of IIK06. The top right panel of Figure [5] shows that the 
density varies by more than two orders of magnitude across the front for all heating rates. 
As \dB/dx\ Xl is increased the density of the warm phase at the far boundary increases and 
hence the temperature decreases. 

For insight into the actual nature of fronts, one may begin by looking at the bulk velocity 
profiles, shown in the bottom left panel of Figure The effect of the initial magnetic field 
strength gradient on the velocity profile of a front is not straightforward. For the lowest 
initial \dB/dx\ Xl case shown the velocity profile is flat and close to zero, as should be the 
case for a static front. As \dB/dx\ Xl is increased the velocity at first becomes larger and 
negative. This is because the saturation pressure is altered from the original hydrodynamic 
value of p sa t/kB = 2612 K cm -3 , as we discussed in §4.1. The models shown here with 
negative velocity profiles are actually condensation fronts. However, as \dB/dx\ Xl is further 
increased there comes a point when the velocity no longer becomes increasingly negative, and 
instead begins to increase. Eventually the front transitions from being a condensation front 
to an evaporation front, which is illustrated by the positive velocity profile of the largest 
\dB/dx\ xl model shown in Figure El This is expected because of the nonmonotonic behavior 
of the saturation pressure as the ambipolar heating rate is increased (see Figure |3]). 

The magnetic field strength profile of the lowest value \dB/dx\ Xl model, given in the 
lower right panel of Figure El is extremely flat. As \dB/dx\ xl is increased the field strength 
decreases across the domain in an almost linear fashion; for sufficiently large values the 
profile becomes nonlinear. Given that the change in density across a front is much more 
dramatic than that of the magnetic field strength, the ratio of the magnetic field strength 
to the number density of neutrals, B/n, changes markedly throughout the transition layer. 

In FigureElwe plot the plasma velocity profile, given by v p ~ v + vp and Equation flBj), of 
each of the front models of Figure [51 The shapes of the profiles are governed by the behavior 
of the magnetic field strength. The plasma velocity is almost constant across the lowest 
\dB/dx\ xl model since the magnetic field strength profile is close to flat, whereas the larger 
\dB/dx\ xl models show more variation in their plasma velocity profiles due to the presence 
of significant gradients in the magnetic field. In all cases the drift velocity is positive and 
larger than the bulk velocity of the front, such that the plasma velocity is also positive. 
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In Figure [7] we compare the ion- neutral drift heating rate, given by Equation (|16p . 
to that of photoelectric heating, given by Equation f|T5|) . for each of the front models of 
Figure O The lowest \dB/dx\ Xl model has a much smaller ambipolar drift heating rate than 
photoelectric heating rate which is why the structure of that front is barely different from the 
hydrodynamic case. The three larger \dB/dx\ Xl fronts have larger ambipolar heating rates 
that are comparable to the photoelectric heating rate. These fronts depart more noticeably 
from the hydrodynamic solution and are less diffusive. 

We also investigate the effect of magnetic field strength on front profiles at fixed initial 
\dB/dx\ Xl . Figure |8] shows a variety of front characteristics for inital field strengths of 1, 3, 
and 5 pG and \dB/dx\ Xl = 308.6 pG pc _1 . The temperature profiles are very similar, with 
the effect of increasing the field strength being larger temperature gradients at small scales 
and thinner fronts. The effect on the density profile is that the higher magnetic field strength 
fronts connect warm phases with higher densities. The velocity profiles are slightly negative, 
which implies that these are actually condensation fronts, and the departure from a static 
solution seems to increase with increasing field strength. Higher field strength solutions have 
flatter magnetic profiles because the efficiency of ambipolar diffusion increases with magnetic 
field strength. 

In Figure |9] we plot the plasma velocity profiles of each of the front models of Figure 
EJ The size of the plasma velocity increases with increasing magnetic field strength, and the 
shape of the profile becomes flatter. This is also due to the higher efficiency of ambipolar 
diffusion at larger magnetic field strengths. 

In Figure [10] we compare the photoelectric and ambipolar heating rates for the front 
models shown in Figure [HJ For these particular cases the photoelectric heating rate is larger 
than the ambipolar heating rate for all field strengths, and the ambipolar heating rate in- 
creases with magnetic field strength. 

4.3. Flux-Freezing Approximation 

For completeness, we also present the flux-freezing approximation, in which the behavior 
of the magnetic field is tied to the density such that B/p is constant in 1-D. This result can 
be obtained by computing the total derivative of the quantity B/p using the continuity and 
induction equations. To calculate the structure of a front in this approximation, we solve 
Equations Q) and ffTUj) and everywhere replace B by pC, where C is a constant. Including the 
definition z = dT/dx, we have a system of three ODEs, which we solve using our shooting 
method, with the mass flux, j, the parameter to be adjusted. We impose the boundary 
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conditions given by equations f|T8|) and f|T9|) . with no need for a condition on the magnetic 
field strength since its behavior is governed by that of the density. 

Figure [IT] shows solutions for an initial density of n = 106.08 cm -3 at various initial 
magnetic field strengths. As the field strength is increased the front becomes thinner and the 
transition reaches a progressively lower temperature, higher density final state at the right- 
hand boundary. Both the density and magnetic field strength span more than two orders of 
magnitude from one phase to the other. While such a range of densities is routinely observed 
in the neutral ISM, such widely varying magnetic field strengths are not (e.g. Troland & 
Heiles 1986), and this provides the first indication that the flux- freezing approximation is 
not suitable for our problem. 

We go on to use these results to calculate ambipolar drift velocities, using Equation (E]), 
and heating rates, given by Equation ffTBT) . throughout the front. These are plotted in the 
lower two panels of Figure [TTJ For the most extreme case shown, a saturated front with an 
initial density of 106.08 cm -3 and an initial magnetic field strength of 5 [iG, we obtain a 
maximum drift velocity of 19.4 km s" 1 and a maximum heating rate of 1.9 x 10 -21 ergs s _1 
cm -3 , three orders of magnitude greater than the photoelectric heating rate. Although the 
equation for drift velocity breaks down for cases in which it is supersonic we may still employ 
it to show that if the flux-freezing approximation held, the drift velocities and heating rates 
would be enormous. Such an outcome is not self-consistent with the rest of the model, and 
allows us to argue that the solutions must be closer to what we have already presented, with 
the magnetic field strength almost constant over the extent of the front for cases in which 
ambipolar drift heating does not dominate. We thus suggest that by the time steady-state 
fronts are established in the neutral ISM the flux-freezing approximation does not apply. 



5. DISCUSSION 

We have shown the magnetic field strength profiles of fronts having ion-neutral drift 
heating rates much smaller than the photoelectric heating rate to be almost flat. In this 
section we argue that it is the thin extent of these fronts that mediates the leakage of 
the magnetic field by ambipolar diffusion. We begin by using our results to justify our 
steady-state and ionization equilibrium assumptions. The minimum flow time through a 
front is of order Tfi ow ~ 0.01 km s _1 /0.1 pc ~ 10 Myr (refer to Figure This should be 
compared to the ion-neutral collision time, given by Tj n ~ (PnTAo) 1 ~ 15.8/n n yr. Thus, 
we have Ti n /rfi ow <C 1 so are safe in our steady-state formulation of ambipolar diffusion. The 
assumption of ionization equilibrium is scrutinized by comparing Tfi ow to the recom bination 



time for hydrogen, given by r rec ~ where ~ 2.06 x 10 U T x l 2 cm 3 s 1 ( Spitzer 
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19781 ). For a front with an initial density of 106.08 cm -3 we calculate a recombination time 
of ~ 70 yr on the cold side, and on the warm side we obtain ~ 5000 yr. For all our other 
front models we also find r rec /rfi ow <C 1, so for this work our simple single-fluid treatment 
of ambipolar diffusion will suffice. 

We now present a diffusive description of fronts in which we compare the thermal and 
ambipolar diffusivities. Taking U = nk B T to be the energy density, we write the thermal 
timescale as r t h = U/pC, and the thermal diffusivity as X t h = kT/U, such that the charac- 
teristic length scale of the problem, the Field length, is given by lp = \J XthJth (Field 1965; 
Begelman & McKee 1990). Hence, the thermal timescale and flow velocity may be written 
in terms of the thermal diffusivity, such that r th ~ lp/Xth an d v t h ~ Kh/^F, respectively. 
In the magnetic field case, the field is redistributed diffusively, with ambipolar diffusivity, 
A^ p = v\r n i, where r n? ; i s the neutral-ion collision time, approximated by r ni ~ 1.58 x 10 3 /rij 



yr (IPadoan et al. Il2000i ). Comparing the thermal and ambipolar diffusivities we obtain: 



A+h K -TljT 1 / 2 , . 

- - 10 ' , 22 



Aad nk B v 2 A r ni B 2 



where is the field strength in units of /iG. We compute Equation (122]) at both boundaries 
of our front models and for all cases obtain tad/toi *C 1. For example, for a front with 
an initial density and magnetic field strength of 106.08 cm -3 and 5 /iG, respectively, and 
\dB/dx\ xl = 308.6 fiG pc _1 , we obtain T AD /r th ~ 4.9 x 10~ 5 on the cold side, and T AD /r th ~ 
4.0 x 10~ 4 on the warm side. This means the drift time is always much smaller than the time 
to flow through the front, suggesting that the field has time to become close to uniform 0. 

Our results show that increasing the ambipolar heating rate changes the structure of 
our front solutions. By balancing the ambipolar and photoelectric heating rates, Equations 
(fT5j) and f|T6|) . and approximating the magnetic field strength gradient as B /L Bcrit , we 
can estimate the critical length scale at which the magnetic field becomes important in 
determining structure: 



J Bcrit 



, 1/2 

Aad -Bo 



47mr PA H 



(23) 



The magnetic length scale is given by L B ~ B/\VB\, so if L B > L Bcrit the effect of the field 
on the structure of a front is small. We compare L B and L Bcr u in Figure [12] for a front with 
an initial density of n ~ 106.08 cm" 3 and initial field strengths of B = 1, 3, and 5 fiG, with 
an initial field strength gradient of \dB/dx\ Xl = 308.6 pG pc" 1 . For the 5 fiG case we obtain 



1 Note that the value of \dB/dx\ Xl enters into this estimate only insofar as it affects the equilibrium 
temperature and front structure. 
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Lb ~ 1.6 x 10~ 2 pc and LbctU ~ 5.4 x 10~ 3 pc on the cold side, and on the warm side we 
find Lb ~ 16.3 pc and Lscrit ~ 2.9 pc. The magnetic length scale is larger than the critical 
scale throughout the front, thus ambipolar drift heating does not have a dramatic effect on 
the structure of a front. 

Previous dynamical studies have claimed that ion-neutral drift is not a sufficiently fast 
diffus ion process for trans porting; magnetic energy, and instead invoked turbulent ambipo lar 
drift (IHeitsch et al. 1120041 ) or turbulent magnetic reconnection ( iSantos-Lima et al. 1120101 ) to 
explain the B — p relation. However, these studies were on larger scales than the fronts con- 
sidered here. Our results suggest that for this simple scenario in which the phase transitions 
are thin, ambipolar diffusion alone is a sufficient mechanism for redistributing the magnetic 
field energy, without the need for turbulence. Our work directly com plements a study o f 



the thermal instability as a formation mechanism for diffuse H I clouds (llnoue et al. 1 120071 ). 



In that work it was shown that ambipolar diffusion is a necessary and sufficient ingredient 
for the formation of a two-phase medium. Once that medium is established the methods 
discussed in this paper may be applied to calculate its structure. 



6. SUMMARY AND CONCLUSIONS 

In this work we have investigated the effect of magnetic fields on two-phase structure in 
the neutral ISM. We have presented a numerical method for calculating the 1-D structure 
of fronts separating the cold neutral medium from the warm neutral medium, including the 
effects of ambipolar diffusion. We showed that the pressure range over which two-phase 
structure is permitted becomes larger, by as much as a factor of two, due to the contribution 
of ambipolar drift heating, with both the minimum and maximum pressures increasing from 
their hydrodynamic values. We find our magnetized front profiles to be very similar to 
the hydrodynamic solutions, and, in cases where photoelectric heating dominates ambipolar 
drift heating, to have close to flat magnetic field strength profiles. We also showed that the 
flux-freezing assumption yields unphysically large drift velocities and frictional heating rates. 
Our method is generic and, by including the appropriate physics, may be extended to other 
astrophysical multi-phase systems. 

Although the 1-D picture discussed in this work is fairly simple, if the magnetic field 
strength and density were related we would have expected to see a correlation. Our results 
are consistent with the observational evidence that there is no relationship between magnetic 
field strength and density in interstellar atomic gas, which suggests that ambipolar diffusion 
is an efficient transport mechanism in the neutral ISM. The effect of ambipolar diffusion on 
the stability properties of thermal fronts will be the subject of forthcoming publications. 
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Fig. 1. — Geometry of front and magnetic field. We seek a front solution connecting the cold 
neutral medium with the warm neutral medium for the case of a uniform magnetic field tangential 
to the front, -E?cnm = B{x)z. The bulk velocity flow is in the x-direction. The warm neutral 
medium quantities to the right of the front are to be solved for. 
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n [cm -3 ] 

Fig. 2. — Thermal equilibrium state of the cooling function, pC(n, T) = 0, in the thermal pressure- 
number density plane. The solid line shows the hydrodynamic case in which there is no ambipolar 
drift heating. The other curves show the effect of increasing the ambipolar drift heating rate at a 
fixed magnetic field strength of 3 pG. In the area above the curves pC > so cooling dominates, 
while below the curves pC < so heating dominates. Increasing the ambipolar drift heating rate 
increases the pressure at which two phases can co-exist. Note that the ambipolar drift heating rate 
is actually a function of density and magnetic field strength, as given by Equation (|16p . but for the 
purposes of this plot we set it to be some constant fraction of the photoelectric heating rate. 
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Fig. 3. — Saturation number density, temperature, and pressure (where Psat/ks = n sa tT sat ) as a 
function of the initial magnetic field strength gradient for initial field strengths of 1, 3, and 5 fiG. 
The range of \dB/dx\ xl has been chosen to show the inflection point of the saturation pressure. 
The 1 //G curve has not been extended further because when \dB/dx\ Xl is too large such a front 
cannot connect to a thermal equilibrium phase before the magnetic field strength becomes negative. 
The 3 and 5 fiG curves can be extended to larger saturation pressures than shown here, until their 
magnetic profiles also become negative. 
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Fig. 4. — Two static fronts with pthermai/^s = 2500 K cm -3 and an initial magnetic field strength 
of 3 juG, but different initial magnetic field strength gradients. Only the phase transition is shown, 
so the reader may find it helpful to picture the cold phase occupying the region x < 0, and the 
warmer medium filling the region beyond the end of the transition. The dashed line shows the 
front subject to a higher ambipolar diffusion heating rate, which connects a cold medium, with 
w-CNM = 84.66 cm" 3 and Tcnm = 29.53 K, with a warmer medium with ji^nm = 0.52 cm" 3 and 
T\vnm = 8063 K. The front with the lower ambipolar drift heating rate connects a cold phase 
having ncNM = 95.51 cm -3 and Tcnm = 25.12 K with a warm phase having ji^nm = 0.32 cm -3 
and Twnm = 8533 K. The top panels show the temperature and density profiles and the lower 
panels show the magnetic field strength profiles and the ratio of the field strength to the density 
for both models. 
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Fig. 5. — Profiles of fronts having an initial density of n = 106.08 cm -3 and an initial magnetic 
field strength of Bq = 3 \xG for various initial magnetic field strength gradients. The top panels 
show the temperature and density profiles and the lower panels show the bulk velocity and magnetic 
field strength profiles of the fronts. The different line styles represent different values of \dB/dx\ Xl 
(note that these same line styles are used in Figures [6] and [7]) . 
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Fig. 6. — Plasma velocity (v p ~ v + vjj) profiles of fronts having an initial density of n = 106.08 
cm~ 3 and an initial magnetic field strength of Bq = 3 fiG for various initial magnetic field strength 
gradients. 
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Fig. 7. — Ambipolar drift heating rates normalized by photoelectric heating rate across fronts 
having an initial density of n = 106.08 cm -3 and an initial magnetic field strength of Bq = 3 fiG, 
for various initial field strength gradients. 
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Fig. 8. — Temperature, density, bulk velocity, and magnetic field strength profiles of fronts having 
an initial density of n = 106.08 cm -3 at various magnetic field strengths for fixed \dB/dx\ xl = 
308.6 fiG pc _1 . The same line styles are also employed in Figures [9l and [TUl 



- 26 - 



0.1 i 1 — i — i i i 1 1 1 1 1 — i — i i 1 1 1 1 1 1 — i — i i 1 1 1 1 1 1 — i — i i 1 1 1 1 1 1 — i — i i 1 1 1 1 




x[pc] 

Fig. 9. — Plasma velocity profiles of fronts having an initial density of n = 106.08 cm -3 at various 
magnetic field strengths for fixed \dB/dx\ xl = 308.6 fiG pc _1 . 
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Fig. 10. — Ambipolar drift heating rates normalized by photoelectric heating rate across fronts 
having an initial density of n = 106.08 cm~ 3 and various initial magnetic field strengths for fixed 
\dB/dx\ Xl = 308.6 fiG pc" 1 . 
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Fig. 11. — Profiles of fronts having an initial density of n = 106.08 cm -3 , hence PthermaiA-B = 2612 
K cm" 3 , calculated in the flux- freezing approximation (without ambipolar drift heating) at various 
magnetic field strengths. The top panels show the temperature and magnetic field strength profiles. 
Although not shown here, the density profile has the same shape as the field strength profile, as 
dictated by flux-freezing. The lower panels show the plasma velocity and the ratio of the ambipolar 
heating rate to the photoelectric heating rate. The solid line shows the hydrodynamic result so 
only appears in the upper left panel. 
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Fig. 12. — Ratio of the critical magnetic length scale to the magnetic length scale for fronts having 
an initial density of n ~ 106.08 cm -3 and various initial magnetic field strengths, with \dB/dx\ xl = 
308.6 /iG pc _1 . Ambipolar drift heating becomes important in determining the structure of the 
front if L B < L Bcr it- 
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Table 1. Properties of cold and warm phases connected by fronts 1 having an initial 
density of tzcnm = 106.08 cm -3 and an initial magnetic field strength of -Bcnm = 3 fiG. 



\dB/dx\ Xl (fiG pc 1 ) 


Front Type 


Tcnm (K) 


"WNM (cm 3 ) 


Twnm (K) 


-Bwnm (a*G) 


Thickness (pc) 


0.31 


Static 


24.63 


0.31 


8580 


3.000 


0.94 


617.2 


Condensation 


25.65 


0.45 


8210 


2.386 


0.28 


1157.3 


Condensation 


27.97 


0.67 


7818 


1.053 


0.12 


1219.0 


Evaporation 


28.30 


0.72 


7743 


0.037 


0.08 



The profiles of the connecting fronts are presented in Figure [5j 



